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ABSTRACT 

The distribution of differential time delays At between images produced by 
strong gravitational lensing contains information on the mass distributions in the 
lensing objects as well as on cosmological parameters such as Ho. We derive an 
explicit expression for the conditional probability distribution function of time 
delays P{At \ 6), given an image separation between multiple images 6, and re¬ 
lated statistics. We consider lensing halos described by the singular isothermal 
sphere (SIS) approximation and by its generalization as proposed by Navarro, 
Frenk, & White (NEW) which has a density profile p oc in the innermost re¬ 
gion. The time delay distribution is very sensitive to these profiles; steeper inner 
slopes tend to produce larger time delays. For example, if Hq = 70 kms“^Mpc“^, 
a A-dominated cosmology and a source redshift cs = 1.27 are assumed, lenses 
with 6 = 5" produce a time delay of Af[yr] = 0.39^q; 22, 0.15^Q;Qg, and 

0.07lj^o;o3g (50% conhdence interval), for SIS, generalized NEW with a = 1.5, 
a = 1.0, and a = 0.5, respectively. At a hxed image separation, the time delay is 
determined by the difference in the lensing potential between the position of the 
two images, which typically occur at different impact parameters. Although the 
values of At are proportional to the inverse of Hq, P{At \ 9) is rather insensitive 
to all other cosmological model parameters, source redshifts, magnification biases 
and so on. A knowledge of P{At \ 6) will also be useful in designing the observ¬ 
ing program of future large scale synoptic variability surveys and for evaluating 
possible selection biases operating against large splitting lens systems. 
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1. Introduction 

The cold dark matter (CDM) scenario predicts relatively cuspy dark matter halos. On 
the basis of systematic N-body simulations, Navarro, Frenk, & White (1996, 1997, here¬ 
after NFW) found that the density prohle obeys the “universal” form p{r) oc r~^{r + rs)“^ 
irrespective of the underlying cosmological parameters, the shape of the primordial fluctua¬ 
tion spectrum and the formation histories. Recent high resolution simulations suggest even 
steeper cusps p oc r~^'^ in the innermost region (Moore et ah 1999; Fukushige & Makino 
2001a,b), and a weak dependence of the inner slope on the halo mass is also reported (Jing 
& Suto 2000). 

The statistics of strong gravitational lensing have been used to probe the density prohles 
of dark halos, e.g., multiple QSO images (Fox & Pen 2001; Keeton & Madau 2001; Wyithe, 
Turner, & Spergel 2001; Li & Ostriker 2001; Takahashi & Chiba 2001) and the long thin 
arcs (Williams, Navarro, & Bartelmann 1999; Meneghetti et al. 2001; Molikawa & Hattori 
2001; Oguri, Taruya, & Suto 2001). These theoretical studies concluded that gravitational 
lensing rates are extremely sensitive to the inner slope of dark halos. The lensing statistics 
of small separations, however, will also be affected by gas cooling and dumpiness in the host 
halo (Keeton 1998; Porciani & Madau 2000; Kochanek & White 2001; Li & Ostriker 2001). 
Thus multiple QSO images with intermediate or large separations, 6^ > S", are more relevant 
in constraining the density prohle of “pure” dark halos. 

At present, several lensing surveys at large separations have been carried out; e.g., the 
Jodrell-Bank VLA Astrometric Survey and the Cosmic Lens All Sky Survey (JVAS/CLASS; 
e.g., Browne & Myers 2001) and Arcminute Radio Cluster-lens Search (ARCS; e.g., Phillips, 
Browne, & Wilkinson 2001). The JVAS/CLASS sample comprises 10,499 radio sources. 
An explicit search for lenses has detected no lenses with image separations Q” < 9 < 15" 
(Phillips et al. 2000), while 18 gravitational lenses with 0.3" < 6 < 3" were found (Helbig 
2000). ARCS also prodnced a null result for lensing events with 15" < 9 < 60" from 1,023 
extended radio sources. The lack of large separation images is only marginally consistent 
if the usnal NFW prohle and A-dominated CDM model are assumed (Li & Ostriker 2001; 
Keeton & Madau 2001), but it also may be ascribed in part to an ehect of the longer time 
delays between more widely separated images. Any intrinsic variability of the QSO will 
result in images which less resemble each other for longer delays (Phillips et al. 2001). 
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So far, time delays between two images have been used primarily to estimate the Hubble 
constant Hq using a detailed lens model of each system (e.g., Grogin & Narayan 1996; 
Barkana et ah 1999). But the importance of mass distribution in estimating Hubble constant 
has been recognized (Impey et ah 1998; Keeton et ah 2000; Witt, Mao, & Keeton 2000) and 
this led to an attempt to constrain the galaxy mass profile from time delays using Monte 
Carlo simulations (Rusin 2000). Instead, in this paper, we consider the statistics of the time 
delay effect analytically. We derive an expression for the cumulative joint probability, i.e., the 
probability that the time delay is larger than At and the image separation is 9, P{>At, 6), 
and also various related statistics. They allow us to estimate the range of probable time 
delays for a given image separation and the extent to which the intrinsic time variability of 
quasars affects observed strong gravitational lensing rates. Our most important result is that 
the distribution of delays is quite sensitive to the density prohles of the lensing objects. For 
example, an NFW density prohle predicts median time delays a factor of three or so smaller 
than the density profile proposed by Moore et al. (1999) and Fukushige & Makino (2001a,b) 
for the same Hq value. While the lenses for which time delays are currently observed are 
dominated by barionic component, we can constrain the density prohle of dark halos if a 
sample of time delays with large separations becomes available. It turns out that the time 
delay statistics are fairly insensitive to other uncertainties, such as the magnihcation bias 
and various cosmological parameters and therefore become a relatively reliable estimator for 
the density prohle of dark halos. 

Of course, all delay values are linearly proportional to the inverse Hubble constant, and 
we here assume its value to be Hq = 70 km s“^Mpc“^, which is consistent with the hnal result 
of Hubble Space Telescope Key Project (Freedman et al. 2001), throughout the remainder 
of our discussion. 

The outline of this paper is as follows. In §2, we briehy describe the usual formulation 
of gravitational lensing statistics. Section 3 presents the analytic formulation of time delay 
statistics. Our main results are shown in §4. Finally we summarize the main results and 
discuss their application in §5. 


2. Description of Gravitational Lensing Statistics 
2.1. Basic Equations 

We denote the image position in the lens plane by ^ and the source position in the 
source plane by ff. We assume spherically symmetric lens objects throughout the paper. In 
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this case, the lens equation (e.g., Schneider, Ehlers, & Falco 1992) reduces to 

y = x-p{x), (1) 

where x = |^|/^o, V = Ihl-^OL/'Co-Dos, ■Co is the characteristic length in the lens plane (see 
§2.3.), and Hol and Dqs denote the angular diameter distances from the observer to the lens 
and the source planes, respectively. The explicit expressions of the scaled deflection angles 
j3{x) for the specihc lens models are presented in §2.3. 

We consider a halo of mass M in the lens plane at a redshift zi^ and a source located 
at ;2s- Then the gravitational lensing cross section a{> 9,> /i), dehned in the lens plane, 
for multiple images with image separation larger than 6 and magnihcation larger than /i is 
given by 

= TT 1/2^0 0(0(M,^s,^l) (2) 

where t/j- is the critical source position to form multiple images (usually given by the position 

of the radial caustic), Q{x) is the Heaviside step function, and p(> p) denotes the fraction 
of area satisfying > /n 

P(>h) = — / dyy<S){p{y)-p). (3) 

y-r Jo 


From this cross section, the probability that a source at ;2s is multiply lensed with image 
separation larger than 9 and magnihcation larger than p is 

/■^s roo 

P{>9,>p-,zs)= dz-L dM7ry‘^^Qp{>n) — {l + zi^)^n{M,zi^), (4) 

Jo JM^in 

where Mmin is determined by solving the equation; 9 = 9{Mrain, zs, zl)- 


To calculate the lens mass distribution, we employ the Press-Schechter function (Press 
& Schechter 1974): 


nps(M,z) 


Po <^o(^) 

VttM < 


d(j M 

1m 


exp 



(5) 


where um is the rms of linear density huctuation on mass scale M at z = 0 and ^ 0 ( 2 :) is the 
critical linear density contrast, (5o(^) ~ 1.69/Zl(z), with D{z) being the linear growth rate 
normalized to unity at 2 : = 0. 


The corresponding differential probability with respect to the image separation is 

d 


P{9,>p-,zs) = 


dzp—{l + zpf 
0 dzp 


dM 

~d£ 




( 6 ) 


M=M„ 
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2.2. Magnification Bias 

So far, we have not considered magnification bias (Turner 1980; Turner, Ostriker, & Gott 
1984). For a source of luminosity L, the effect of the magnihcation bias may be included as 
follows: 

1 /‘°° d 1 

P^{0]Zs,L) = — -— / dn —P{e,>fi]zs) ^zs,L/fi)- 

^{zs,L)Ji d/i /i 

, Cdt,^ ,o fdM 9 .O r\ /^,\ 

= / d2;L;^(l + ^L) —TTy^ilB{zs,L)n{M,zi^) , (7) 

j 0 L J M=Mmin 

where <F( 2 :s,T) is the luminosity function of sources and B{zs,L) is 

Thus the probability distribution taking account of the magnihcation bias is simply expressed 
as equation (6) with p(> fj,) replaced by B{zs,L). By neglecting magnihcation bias for 
simplicity, we may simply set B{zs,L) = 1. 

It should be noted that fi{y) may be interpreted as the magnihcation of the total im¬ 
ages, of the brighter, or of the fainter image among the outer two images, depending on the 
observational selection procedure for hnding gravitational lens systems (Sasaki & Takahara 
1993; Cen et al. 1994). As indicated in §2.3.2, however, such a diherent choice of magnih¬ 
cation does not ahect the conditional time delay probability very much, especially for small 
inner slope values, while it changes the lensing rate by a factor of two. Thus we choose to 
designate the magnihcation fi{y) in the bias factor (eq. [8]) as the total magnihcation of all 
images throughout this paper. 


2.3. Specific Density Profiles 
2.3.1. Singular Isothermal Sphere 
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Then the lens equation has two solutions x± = y \i \y\ <yr = 1- The separation between 
two image may be written as 


^ ^ 0 ( 3 ;+ - x_) 
-Dol 

The magnihcation of each image is 




y dy 
x± dx± 


and their total magnihcation is given by 




Stt 


^LS 

c/ Bos' 



1 

-± 1 , 

y 


2 

y' 


( 11 ) 


( 12 ) 


(13) 


To compute the probability distribution functions ( 6 ) and (7), we also convert the mass 
function (eq. [5]) to a velocity function by using the spherical collapse model (e.g., Nakamura 
& Suto 1997). 


2.3.2. Generalized NFW Profile 


The halo density prohles predicted by recent N-body simulations may be parameterized 
as a one-parameter family, the generalized NFW prohle (Jing & Suto 2000): 


p(r) 


Pcrit<^c 

(r/rg)“ (1 -hr/rs)^““’ 


(14) 


where r^ is a scale radius and 5c is a characteristic density. The prohle with a = 1 corresponds 
that NFW originally proposed, while the prohle with a = 1.5 resembles the one claimed by 
Moore et al. (1999) and Fukushige & Makino (2001a,b). The scale radius generally depends 
on the mass and the redshift, and is related to the concentration parameter: 


Cvir(Tf, z'^ 

The characteristic density is given by 


r^ifiM,z) 
rfiM,z) ■ 


Ayii-flvir _ (3 - a)c"^. _ 

3 2 F 1 (3 — a, 3 — a; 4 — a; —cvir) ’ 


(15) 


(16) 


with 2 F 1 (a, b] c; x) being the hypergeometric function (e.g., Keeton & Madau 2001). Follow¬ 
ing Oguri et al. (2001), we adopt the mass and redshift dependence reported by Bullock et 
al. ( 2001 ) and consider the median amplitude of the concentration parameter Cnorm: 


Cvir(Tf, z'j Cn 


2 — a 


M 


1 + 2: V lO^^h'^M© 


- 0.13 


( 17 ) 
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where h denotes the Hubble constant in units of 100knis“^Mpc“^. Scatter of the concentra¬ 
tion parameter is modeled by a log-normal function with the dispersion of cTc = 0.18 (Jing 
2000; Bullock et ah 2001). We £x Cnorm to the value estimated in the simulations (Bullock 
et ah 2001), Cnorm = 8, throughout this paper. 


In these profiles, we choose .^o = ^'s and denote the scaled deflection angle I3{x) as 
I3{x) = b f{x). The dimensionless factor b (of order unity) and the function f{x) are related 
to the dark halo profile as follows: 


fix) 


b = 


lOvrGpcrit^crs DqlDls 


D, 


os 


1 

X 



X 


^1 -I- 


3 — 0 ; * 


(18) 

(19) 


The lens equation has three solutions xi,X 2 ,X 3 {xi > X 2 > x^) if \y\ < i/r, where Ur is 
the position of the radial caustic. The image separation is defined between the outer two 
solution and is approximated as (Hinshaw & Krauss 1987): 


. _ ^o(a;i - xs) 2r^xt 
u = - ~ -, 

-Dol Doh 


( 20 ) 


where Xt is the position of the tangential critical curve (the Einstein radius). The approx¬ 
imation (eq. [20]) is sufficiently accurate for the range of interest here (Upper panels of 
Figure 1), and is useful in computing the lensing probability. 


The top panels of Figure 2 shows that the total magnification for the generalized NFW 
profile is well approximated as (Blandford & Narayan 1986; Kovner 1987; Nakamura 1996) 


y{y) ~ 


Vr 

htO 

y 


hrO 


(1 -y/l/r)^/^ 


iy ^ l/crit)) 

iy > l/crit); 


( 21 ) 


where pto and pro are 


ytQ — 


^.Xt 

yril - bf'{xt)y 


Xr 


hrO 1 


yr V yrbf'ix,) ’ 

with x-r being the position of the radial critical curve. Finally |/crit is given by 

!/crlt _ + Vi'm + 411,V? 


( 22 ) 

(23) 


,2 

^rO 


l/r 




(24) 



The magnification of the brighter or fainter image can be approximated as follows: 


P'bright(|/) — 


htOl/r 

2 y 

hrO 1 

2 {l-y/yrY/"^ 


{y ^ l/crit)) 


{y ^ l/crit)) 


(25) 


P'faint(|/) — 


yto yr 
2 !/■ 


(26) 


Figure 2 shows the comparison between the magnihcations of the three cases and their 
approximation described above (eqs. [21], [25], and [26]). Figure 2 indicates that the above 
approximation breaks down around its minimum value. In practice, however, this level of 
discrepancy does not affect the result of the magnihcation bias (Nakamura 1996). More 
importantly, the difference in magnihcation using total, brighter, or fainter images is only 
about a factor of two, as easily seen from the above approximations. That is, the different 
choices of the magnihcation only increase (or decrease) the probability by a factor of two 
for all separations. This numerical factor is canceled out in calculating the conditional 
probability of time delays because it is dehned by the ratio of these probabilities (see §3.1). In 
the SIS case, however, the diherence between the magnihcation of brighter and fainter image 
is not negligible (see eq. [12]) and as a result the ehect of diherent choice of magnihcation on 
time delay probability distributions becomes larger. Therefore, we conclude that the ehect 
of diherent dehnitions of the magnihcation is small in time delay statistics, especially for a 
small inner slope {a < 1.5). 


3. Analytic Formulation of the Time Delay Probability Distribution 

3.1. General Formulation 

The lens alters the time taken for light to reach the observers and inevitably produces a 
diherential time delay between multiple images (Refsdal 1964, 1966). The value of the time 
delay is calculated as (e.g., Schneider et al. 1992) 

cAt(|/) = -^2^^(1 + ^l) [(l){x^^\y) - (l){x^‘^\y)] , (27) 

-^OL-^LS 

where x^^'^ {i = 1,2) are two image positions. The Fermat potential is written as 

Y{x,y) = (28) 

in terms of the lensing potential The explicit expressions for the time delays in the 

SIS and the generalized NFW prohles are given in §3.2. 
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The cumulative probability of time delays is calculated from equations (7), (8), and 


(27): 


P^{>At,e;zs,L)= I dz-L -^{1 + z-tf 


'0 


dM 

ltd 


TT l/r ^0 ^s, T)n(M, ^l) 


M=M„ 


{> At; zs, L) = { y^^{zs,L) J 

J/min 

0 


J/r 

dyyHzs,L/fi{y))—- {y 

min l/r)? 




(ymin ^ yr)i 


( 29 ) 

( 30 ) 


where the superscript T on the magnihcation factor means that the time delay threshold is 
taken into account. The lower limit of integral i/min is determined by solving the equation; 
At = At(|/min). This expression is valid when the time delay is a monotonic function of the 
source position y as in all our examples below. 

When the magnihcation bias is neglected, we can replace B^{> At; 2 : 3 , L) by 




yn 


Afy>A() = <!y?/ ^9!' 1 (fcinO.). (31) 


0 


(jjmin ^ l/r 


The cumulative conditional probability can be calculated from equations (6) and (29): 

P^{>At,e;zs,L) 


Pl>At\e;zs,L) = 


P^{0;zs,L) 


(32) 


We dehne the joint probability distribution of time delays as: 


P^{At,9;zs,L) = 


d 


^(^^-p=(>AM;.s.L) 

/ dzi ,-—(1 + 2:l) 

Iq dz-L 


^ B'^{At; Zs, L)n{M, zl) 


B^{At;zs,L) = < y^^{zs,L) 


y^{zs,L/y{y)) 


d{At) 


y{y) 


J/=J/min 


0 



(33) 

J 

< yv), 

(34) 

> yr)- 



Similarly, when the magnihcation bias is not taken into account, we replace B'^{At; zs, L) 
by 


f ^ 

dy 

N^{At) = { y? 

_d{Atf _ 

lo 



y=yinin 


(l/min < l/r)) 
(jjvciin ^ l/r)- 


(35) 
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Finally the conditional probability distribution of time delays is calculated as in equation 


(32): 




P^(At,0;zs,L) 

P®(e;zs,i) 


(36) 


3.2. SIS and Generalized NFW Profiles 

The time delay in the SIS case is calculated from the lensing potential, 'ip{x) = |a;|, and 
equation (27): 

cAt{y) = 3271^ f-] ^ 

Vc/ Dos 

In the generalized NFW case, it is useful to adopt the following approximation in cal¬ 
culating the time delay: 

(/>(x3,y) - ^ 2xty. (38) 

The accuracy of this approximation is shown in Figure 1 {Lower panels). In the analysis 
below, we employ the approximations shown in Figures 1 and 2. From equations (27) and 
(38), the time delay for the generalized NFW case is given by 

c^tiy) = + ^l)|/- (39) 

-^OL-^LS 


4. Results 

4.1. Setting the Conditions 

In what follows, we consider three representative cosmological models dominated by 
CDM; Lambda CDM (LCDM) with (Hq;^W s) = (0.3,0.7,0.7,1.04), Standard CDM 
(SCDM) with (1.0,0.0,0.7,0.54), and Open CDM (OCDM) with (0.45,0.0,0.7,0.83). The 
amplitude of the mass fluctuation, as, is normalized so as to reproduce the X-ray luminosity 
and temperature functions of clusters (Kitayama & Suto 1997). 

We calculate the probability distribution of image separations and time delays adopting 
the detection condition in the JVAS/CLASS survey. The sample of radio sources have a flux 
distribution with ‘h(S') oc (Rusin & Tegmark 2001). In this case, magnification bias 

factors (8) and (30) are simplihed and can be written as 

dyy {y{y)y~\ 



B{zs,L) 


(40) 



and 
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2 

— dyy (l/min < l/r), 

dr JymATL 

0 (l/min ^ Hr) ■ 


(41) 


with 7 being the slope of the radio Inminosity fnnction, 7 = 2.1. Using the fact that the 
Inminosity fnnction is described by a power-law and the time delay is proportional to the 
sonrce position (i.e., At oc y), equations (34) and (35) are also simplihed as 


and 


B^{At; zs,L) 


Aty^ 

0 (l/min ^ yr)i 


N^{At) 


2|/min 

At 1/2 
0 


{Vrain ^ l/r); 
{Vrain ^ Ur)- 


(42) 


(43) 


The redshift distribution of the parent population of radio sources is not known in the 
JVAS/CLASS survey. Thus we £x the source redshift to the mean redshift of a 27 object 
subsample, (zq) = 1.27 (Marlow et ah 2000). The results are not substantially different if 
one takes into account the observed redshift distribution of the subsample (Keeton & Madau 
2001 ). 


4.2. Time Delay Probability Distribution 

Figure 3 plots the predicted probability distribution of image separations for various 
density profiles {Left panel)] SIS, generalized NFW with a = 1.5, a = 1.0, and a = 0.5, and 
also for various cosmological models {Right panel) fixing a = 1.5, for reference. As shown in 
these plots, the probability strongly depends on the inner profile of dark halos. The steeper 
inner profile produces multiple lenses much more efficiently. The JVAS/CLASS data are 
also shown in Figure 3. This plot indicates that small separation lensing is consistent with 
the SIS model, while the SIS model predicts too much large separation lensing, as previously 
shown by Li & Ostriker (2001). 

Consider next the time delay probability distributions. Figure 4 plots the cumulative 
conditional probability (eq. [32]) for different density profiles, and Figure 5 plots its loga¬ 
rithmic time derivative: 

P{> At I 6] zs, L) = P{At I 9] zs, L)At. 


d{\nAt) 


(44) 
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The cosmological model is LCDM in each case. We present results both with and without 
the magnihcation bias effects. In the Top panel of Figure 4, we plot the three cases of 
magnihcation bias dehned by that of the total, brighter, and fainter images. The effect of 
different choice of magnihcation in the generalized NFW case is smaller than in the SIS 
case, as discussed in §2.3.2. These plots indicate that time delay statistics depend strongly 
on the inner slope of the density prohle, a. Steeper inner prohles tend to produce longer 
time delays. This is probably because of the deeper gravitational potential for steeper inner 
prohles. The non-geometric part of the time delay is simply proportional to the diherence 
in the potential value at the image positions. Therefore, we can constrain the density prohle 
of dark halos from this statistics if a reliable and unbiased sample of time delays becomes 
available. 

We note the diherence in the asymptotic behavior in the limit of small At for diherent 
magnihcation biases. For the SIS prohle, we obtain the asymptotic behavior of the condi¬ 
tional probability distribution in the limit of the small At from equations (13), (37) and 
(42): 

P{At\e;zs,L) (X {Atf-'^. (45) 

For the generalized NFW prohle, we also obtain the same relation from equations (21), (39) 
and (42). The asymptotic behavior (45) indicates that the uncertainty in the slope of the 
QSO luminosity function severely ahects the probability of shorter time delays. We are 
not concerned with details of P{At \ 9) at small At, where various uncertainties such as the 
hnite size of the sources may ahect the shape of P{At \ 9) anyway. Except for the small 
At regime, however, the uncertainty of the magnihcation bias is negligible and the time 
delay probability distribution depends primarily on the inner slope of lens density prohle. In 
fact, the cumulative conditional probability P{> At\9) with or without magnihcation bias 
is quite similar and probably observationally indistinguishable for the foreseeable future. 
In particular, the insensitivity of magnihcation bias shown in Figures 4 and 5 is in marked 
contrast with the usual lensing probability as a function of image separation which is ahected 
by more than one order of magnitude (Wyithe et al. 2001; Li & Ostriker 2001; Takahashi & 
Chiba 2001). 

To see the dependence of separations and density prohles on the conditional probability 
distributions more clearly, we calculate the median time delays At^ed as follows: 

P(>Atmed I 6';2 :s,T) = ^. (46) 

We plot the Aimed as a function of separations 9 in Figure 6. The comparison with the 
observation in this hgure will be discussed in §4.3. The error-bars in this plot are dehned 
by the ±25% level. This hgure shows that the median time delays are well htted by the 
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power-law: 

At^ed = (47) 

The best-fit values for A and B are summarized in Table 1. Since ~ 2 for all density 
profiles, Figure 7 plots Atmed/^^ to illustrate their strong dependence on a. It is clear in 
this hgure that the difference of Aimed between prohles with various a values is larger for 
smaller separations. To probe the density prohle of dark halos, however, small separation 
lenses may suffer from the complex physics on baryonic components. Therefore lenses with 
S" ^ ^ ^ lO", which are typically associated with galaxy gronps or clusters, may be more 
useful in constraining the density prohle of dark halos. Figure 7 also indicates that it is 
difhcnlt to determine the inner slope from a single observed lensing system, and several 
leasing systems are required to reduce the statistical uncertainty. 

Next consider the dependence of the statistics on cosmological models or sonrce redshift. 
Figure 8 displays the predictions for different models and sonrce redshifts, indicating that 
their effects are rather small. In particular, we emphasize that the time delay probability 
distribution is insensitive to ag, while the SCDM model with ag = 1 yields nearly two orders 
more mnltiple lenses compared with ug = 0.54. The time delay probability distribntion is 
also insensitive to the redshift uncertainty, especially at high redshift ; 2 g > 1. Therefore, 
uncertainties of both cosmological models and redshift are not a signihcant source of uncer¬ 
tainty in time delay statistics. This robustness is also an advantage of time delay statistics 
compared with the usnal overall lensing rate statistic. 


4.3. Comparison with Existing Observational Data 

Here we tentatively compare our theoretical results with the existing observational data, 
although the detailed modeling is already available for each system in our small sample. The 
adopted data are summarized in Table 2. We use time delay data for two images systems 
and exclude the four image lens systems for simplicity. The data points in Figure 6 and 
Figure 9 show the comparison with observations for fonr lens systems. Although Figure 9 
displays a more direct comparison. Figure 6 may be useful also in examining whether there 
is any dependence of the inner slope a on the mass of the lensing objects. For Figure 6, we 
hx the source redshift at 2:3 = 1.27, although the source redshift for each system is already 
known. This does not make any substantial differences as indicated in Figure 8. We plot 
the theoretical predictions of SIS and generalized NFW density prohles. Observational time 
delays are indicated by triangles (in Fig. 6) and arrows (in Fig. 9). From these plots, we 
hnd that the existing data prefer the steeper inner prohles, SIS prohle or generalized NFW 
prohle with a = 1.5. Although a diherent choice of magnihcation ahects the median time 
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delays by at most a factor of two (see Fig. 4) in the SIS case and less in the generalized NFW 
case, this is smaller than the error-bars of median time delays and does not change our basic 
conclusion. The steep density prohles we infer are also consistent with the detailed models 
(e.g., Keeton et al. 2000) of these particular lens systems. Of course, the image separations 
in the lens sample considered here correspond to the inner regions of the lensing galaxies 
where a plethora of evidence (especially rotation curves) already lead us to expect a roughly 
isothermal distribution of the total mass, luminous and dark combined, so the results are 
hardly surprising. 


5. Summary and Discussions 

In this paper, we formulated the differential time delay probability distribution and ex¬ 
amined its dependence on lens density prohles as well as on magnihcation bias, cosmological 
model, and source redshift. We found that the probability distribution of time delays de¬ 
pends most sensitively on the inner slope of the lens density prohles. The diherence between 
the various density prohles we examined is very large, more than one order. For example, if 
Hq = 70km/s/Mpc, A-dominated cosmology and a source redshift 2:3 = 1.27 are assumed, 
lenses with 6 = 5" induce the time delay; At[yr] = 1.5lg;g, 0 . 39 lo; 22 , 0.15lg;o9, and 0.07llg;o|g 
(50% level), for SIS, generalized NFW with a = 1.5, a = 1.0, and a = 0.5, respectively. 

On the other hand, the ehects of cosmological models, source redshifts or magnihcation 
bias are rather small, aside from the well known and important linear dependence on the 
inverse Hubble constant. Moreover, the delay distributions are quite insensitive to the nor¬ 
malization of the overall lensing rate, because the conditional probability is dehned by the 
joint probability divided by the usual lensing probability and the dihering of normalizations 
are almost canceled out. Therefore, one could strongly constrain the core structure of dark 
halos if a large sample of time delays becomes available in future systematic survey. 

Although existing lens systems are individually modeled in detail including the central 
density prohles, such careful treatment of each lens system in upcoming very large samples 
may not always be practical. Thus a statistical treatment of time delays is informative. For 
example, we can predict the range of probable time delays of large separation lenses from this 
formulation if the density prohle of dark halos is fully settled. This would allow estimation 
of a plausible range of time delays for future lens systems even when the lensing object has 
not been identihed. 

Comparison with the meager existing observational data suggests that the density prohle 
has a rather steep cusp, a ~ 2 , although it is already known that the current observed systems 
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are dominated by baryonic component. Moreover, it seems that the smaller separation lenses 
prefer steeper inner prohles. This result may support a two-population model, i.e., galactic 
mass halos with a steep inner slope {a ~ 2) and cluster mass halos with a shallower inner 
slope (a < 1.5), as proposed by Keeton (1998) (see also Porciani & Madau 2000; Li & 
Ostriker 2001). To see whether there are two (or more) populations of halos, time delay 
data for various separations, especially separations with 5" < 6* < lO", will be essential. 

Another application of time delay statistics is that future large samples of lenses will 
probably be systematically monitored by one or more of the ambitious synoptic surveys now 
being planned. Designing an efficient sampling rate and observing strategy requires an idea 
of the range of time delays that might reasonably be expected. 
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Table 1. Fitting Parameters of Median Time Delays 


Model 

71 [yr] 

B 

SIS 

8.05 X 10-2 

1.77 

q;=1.5 

1.91 X 10-2 

1.86 

q;=1.0 

5.37 X 10-3 

2.05 

a=0.5 

1.48 X 10-3 

2.38 


Note. — Fitting Parame¬ 
ters A and B are defined in 
equation (47) 


Table 2. Observed Time Delays for Lenses with Two Images 


Lens 

2:s 

At [yr] 

0 

Ref. 

B0218+357 

0.96 

0.028 ± 0.004 

0.33" 

1,2 

PKS1830-211 

2.51 

0.071 ± 0.014 

0.97" 

3,4 

B1600+434 

1.59 

0.139 ± 0.011 

1.39" 

5, 6 

Q0957+561 

1.41 

1.157 ± 0.002 

6.17" 

7,8 


References. — (1) Cohen et al. 2000; (2) Biggs et 
ah 1999; (3) Lovell et al. 1998; (4) Lehar et al. 2000; 
(5) Burud et al. 2000; (6) Koopmans et al. 2000; (7) 
Oscoz et al. 2001; (8) Kundic et al. 1997 
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Fig. 1.— Accuracy of approximations for generalized NFW profiles with various parameters 
b (eq. [18]) and the inner slope a] a = 1.0 {left panels) and a = 1.5 {right panels). Upper 
panels display the ratio of the image separation and its approximation (eq. [20]). Lower 
panels give the ratio of the difference of the Fermat potential and its approximation (eq. 
138]). 
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Fig. 2.— The magnification and its approximation defined in §2.3.2 for generalized NFW 
profiles with various parameters b and the inner slope a. From top to bottom, the magnih- 
cation is dehned by that of the total images, of the brighter and of the fainter image among 
the outer two images, respectively. 
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0[arcsec] 


Fig. 3.— Differential Probability of lensing in the JVAS/CLASS survey. The magnihcation 
bias is included. Left panel shows the probability for various density prohles; SIS {solid), 
generalized NFW with a = 1.5 {dotted), a = 1.0 {short dashed), and a = 0.5 {short dash- 
dotted). Right panelis the plot for various cosmological models; LCDM {dotted), SCDM {long 
dashed), and OCDM {long dash-dotted). The observational value is shown by the histogram. 
The la constraint from the JVAS/CLASS null result is also shown by a horizontal line with 
a downward arrow. 
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Fig. 4.— The cumulative conditional probability of time delays (eq. [32]) for various density 
profiles with image separations 6 = l" (solid), 5" (dotted), lO" (dashed), and 30" (dash- 
dotted). The cosmological model is LCDM in each case. Results, with (thick) and without 
(thin) the magnihcation bias, are presented. At the top panel, the effect of different choice of 
magnihcation is shown by the same three lines and shadings; magnihcation is dehned by that 
of total images (center lines), of the brighter (right lines) and of the fainter image among 
the outer two images (left lines). 
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Fig. 5.— Same as Figure 4, but the differential conditional probability distribution of time 
delays times At (eq. [44]) is plotted. 
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Fig. 6.— Median time delays, Aimed, as a function of separations 9, for various density 
profiles. The median is calculated by equation (46). The error-bars in this plot are dehned 
by the ±25% level. A line with At oc 6^ is also shown for reference. Some observational data 
(see Table 2) are indicated by triangles. 
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Fig. 7.— Dependence of time delays on inner slope a. Median time delays divided by the 
square of the separations, Atmed/6*^, are plotted. In this plot, the SIS is regarded as a = 2.0. 
The error-bars are the same as in Figure 6. 
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Fig. 8.— Dependence of the conditional probability distribntion of time delays on cos¬ 
mological model and sonrce redshift. A generalized NFW density profile with a = 1.5 is 
adopted. In the upper panel, the cosmological models are LCDM (solid) , SCDM (dotted), 
and OCDM (dashed). Dash-dotted lines are also for SCDM case, but in this case erg = 1 
instead of cjg = 0.54. In the lower panel, source redshifts zs = 0.5 (solid), zs = 1.0 (dotted), 
zs = 1-5 (dashed), and zg = 2.0 (dash-dotted) are plotted. 
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Fig. 9.— Comparison of predictions with existing time delay data using a cumulative con¬ 
ditional probability. The sample of the lens systems used here is summarized in Table 2. 
Arrows show the observed values of the time delay in each system. Theoretical predic¬ 
tions with various density prohles are displayed by lines; SIS {solid), generalized NFW with 
a = 1.5 {dotted), a = 1.0 {dashed), and a = 0.5 {dash-dotted). In calculating the theoretical 
distributions, the LCDM model was assumed. 


















